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The parton model for semileptonic B meson decays is studied with 
special attention to the decay distributions. We find that the spectra 
show dramatic variations when we introduce cuts on the hadronic 
energy or invariant mass of hadrons. Results for both 6 — > n and 6 — > c 
decays are presented. The detailed spectra may help to separate the 
two types of decays. 
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The decays of B-mesons have been studied extensively and have been very 
useful in extracting properties of the weak interactions. In particular, the 
b ^ u transition is still an active field of investigation. On the experimental 
side there is a large effort hindered by the difficulty of identifying the decays 
of B-mesons to light quarks. This difficulty of identifying h ^ u events forced 
investigators to concentrate on the end-point spectrum of the electron energy 
or study exclusive decaysfll]. In this letter we present an extended study of 
the semileptonic B-decays. We are interested in the decay 

B ^ Xu + e-y (1) 

for which we wish to give explicit distributions in several kinematic variables. 
With the decay spectrum completely determined by a single parameter in 
the distribution function, it may be compared to the competing process in 
B — i> Xc + e~z/. With the help of these distributions it should be possible 
to test characteristics of the decay products with various kinematical cuts, 
and, in particular, collect events typical for this process. The hope is that 
with the results presented here the experimental analyses could incorporate 
many more events than those confined to the endpoint region of the electron 
spectrum. 

There are three models for the inclusive decays of reaction (1). The early 
model of Altarelli et al. [@], to be denoted as ACM, uses a spectator model and 
a distribution of quarks within the B-meson described by the Fermi motion of 
the spectator quark. This approach treats the phase space effects correctly, 
but is rather crude, as the authors state^, because it depends on an unknown 
distribution function for the spectator quark. The model pictures the overall 
decay as the disintegration of the B-meson into the spectator quark plus the 
decay products of the heavy quark. 

The second approach]^, ^ uses the parton model in an infinite momentum 
frame (IMF). The probability of finding a b-quark in a B-meson carrying a 
fraction x of the mesons momentum in the IMF is given by the distribution 
function f{x). The distribution function f{x) has the functional form sug- 
gested by theoretical arguments]^, H, 0] and it peaks at large values of x. 
The kinematics for the decay of a heavy to a light quark involve the corre- 
lation of two currents at short distances, where the incoherence of the decay 
products is justified. This model has recently been improved and applied 



to double differential inclusive decay distributions in reference Q where a 
detailed treatment of the kinematics and the allowed physical regions can be 
found. A special feature of the parton model is the distribution of b quark 
in the B-meson, which provides a continuum spectrum for the mass of the 
recoiling hadrons. Here we only summarize the main results and point out 
new features before turning to applications. 

A third model, suggested very recently 0, improves the free quark model 
by including I/ttlq corrections from the heavy quark theory. These correc- 
tions turn out to be rather small. Obviously, the three models incorporate 
dynamical corrections in a different manner. In the parton model dynamic 
corrections are included through the experimentally determined distribution 
function and contributions from two scaling variables x+ and x_, to be de- 
scribed below. 

In the parton model 0, ^ the decay kinematics are different from those 
in deep inelastic scattering. As a consequence, the energy-momentum con- 
serving (5-function has now two roots, that is, it gives two scaling variables 

which are the light-cone variables of the current-momentum. Thus we have 
again scaling of the distribution function but now in terms of two variables. 
We define the kinematic variables: Pb= momentum of B-meson; Pe,Ee = 
momentum, energy of electron; P^, = momentum of neutrino; Px = momen- 
tum of hadrons; q = Pe + P^ = momentum of current; and Mx = invariant 
mass of the final hadronic system. In terms of these the decay is in general 
defined as: 

^^ ~ {2nfMB "'2E, 2E/ ^^^ 

The leptonic tensor has the simple form 

The hadronic tensor is defined in analogy to deep inelastic scattering, 
W^,{PB,q) = -g,uW,{q^q■PB)+^^^^W2{q^q■PB)-^e^uaf3^Ws{q^q■PB)+■ 

B B 

(5) 



where the integration over the hadronic variables has aheady been performed. 
After carrying out the integrations over x we arrive at the structure functions^ 

W,{q\q-PB) = 2[/(x+)-/(x_)], (6) 

W2iq',q-PB)/Ml = ^^[x+/(a;+)+x_/(x_)], (7) 

Ws{q',q-PB)/Ml = -^-^[/(x+) + /(x_)]. (8) 

Mb \q\ 

and hence the h —^ u triple differential decay rate 
dV G''\Vui\\^-E, 



dEedq^ds Stt^Mb \q\ 



{x+f{x+){2Ee - MbxJ) + (x+ ^ x_)}. (9) 



This formula is simple and shows the dependence on the two light-cone vari- 
ables x± . 

The distribution function f{x) can be taken from references |0, P, |3l- Its 
functional form is very similar to the fragmentation function of a b quark into 
a B-meson. It has also been argued, on physical grounds, that the distribu- 
tion and fragmentation functions for heavy mesons are the same (reciprocity 
relation) p, |10], |Tl|. The fragmentation function was measured in several 



experiments and can be represented [jTT], |T2| as 

with Ep a parameter and N the normalization factor. We will use three values 
£p = 0.003, 0.006. 

Finally, we must comment on the kinematic regions where the model is 
valid. Two criteria must be fulfilled. 

(i) The decay involves the correlation of two currents. We are more 
confident in applying the model when this distance between the two currents 
is close to the light-cone or at short distances. 



^Note that the structure functions in the present article are different from those in 



(ii) The model should apply in the region where many particles are pro- 
duced. Thus, we must avoid the edge of the phase space where one or two 
pions are produced. 

To sum up, we feel rather confident in using the parton model for the 
decay B -^ X^ + e~z/ away from the edges of the phase space. For the decay 
B -^ Xc + e^u, we are less confident, because the distances involved are 
larger. Furthermore, for the decays B ^ Xc + e~u we must keep the c-quark 
mass, which modifies equation (9), and it makes a difference whether we 
choose the minimum value of the final hadronic mass as rric or Md, as we 
will show in figures lb and 2a. 

As we mentioned, one of the purposes of this work is to present detailed 
spectra which can eventually be compared with experiments. The simplest 
parameter to measure is the electron energy. In addition one may be able to 
measure the total hadronic energy Ex = Mb — qo- In Fig. la we show the 
h ^ u double differential decay rate in E^. and go- The spectra show a striking 
dependence of go- Most of the events occur for O.SM^ < Ex < CSM^. 
For smaller hadronic energies the spectra shift to higher electron energies, 
Ee > 2.0 GeV . This correlation of events may help to isolate h ^ u events. 

For comparison we also calculated the distribution for the B ^ Xc + e^v 
decay. In Fig. lb we show the double differential decay rate in E^ and go. 
Comparing with Fig. la we note that this channel runs out of events at 
Ef. = 2.25 GeV and the events above 2 GeV are very few. By making a cut 
in go > Mb — Mb = 0.65Mb there are no events left for b -^ c. This can be 
used as another criterion for isolating b ^ u events. 

The integrated spectrum dV/dEe is shown in Fig. 2a for two fragmentation 
parameters Sp. In the same figure we show the corresponding spectra for 
b ^ c. For these curves we have chosen \Vch\ = \Vub\- It is clear that the 
b ^ c spectrum is softer and a very small fraction lies at E^. > 2.0 GeV. We 
show two spectra for the b ^ c decays, by varying the minimum value of the 
final hadronic mass. We have chosen two values {Mx)min = f^c = 1.5 GeV 
and {Mx)min = Mb = 1.86 GeV. We notice that the total decay rate varies 
considerably when we vary {Mx)min- We feel that the final quark mass should 



be replaced by the running charm quark massQ. As a result, the choice of 
{Mx)min influences the determinations of \Vcb\- 

For comparison we calculated the b ^ u spectrum dT/dEg in the ACM 
model0]. We show in Fig. 2b the ACM and parton spectra together. (In 
these spectra we include QCD radiative corrections which appear as a multi- 
plicative factor, as described in p after the work of p, |13[.) We note that the 
parton spectrum is lower and falls off less steeply than in the ACM model. 
Consequently, the values for Vub extracted from the end-point energy will be 
larger in the parton model with the experimental fragmentation function of 
eqn (10). Of course we can reproduce the free quark model by modifying 
the fragmentation function so that it has the limit 5(1 — x) when the width 
parameter tends to zero, but this is not what is being measured for the frag- 
mentation function |[Tl|, |12|, 0. This means that the semileptonic rate of a 
free b quark is substantially greater than that of a 6 quark confined in a -B 
meson. This is analogous to deep inelastic scattering, where moments of the 



distribution functions were measured to be smaller than one 14 . Thus the 



spectra in Fig. 2b differ not only in shape but also in overall scale. 

The double differential b —>■ u decay rate dT/dEgdMx depends strongly 
on the invariant mass of the hadrons Mx- We show in Fig. 3 dT/dEedMx for 
various values of Mx- The decay spectrum is larger in the range 0.2Mb < 
Mx < OAMb- The characteristic feature again is that the spectra shift to 
larger values of E^. as Mx decreases. After integration over Eg the sum of 
the curves gives the mass distribution, which agrees with mass distributions 
published before [Q and is similar but not identical with the curves in ref. 



The parton model provides an interesting alternative for analysing inclu- 
sive semileptonic decays in important kinematical regions where the decay 
is dominated by short distance physics. We have produced decay spectra 
for a variety of interesting and physically accessible observables which probe 
the decay dynamics in a much more complete way than the simple electron 
spectrum. These spectra are obtained from a simple and compact formula 
(eq. (9)) based on a one parameter parton model. As high statistics data 



*This topic is still under investigation. 



become available this model may play an increasingly important role in sep- 
arating h ^ u from b ^ c decays and lead to a better determination of the 
Vub/Vcb- In fact, the distribution function f{x) can be measured, in principle, 
in the b ^ u decays and then be used as an input to calculate the b ^ c 
spectra. 
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Figure Captions 



la. The b ^ u double differential decay rate dT / dEedqo vs Ee for various 
values of go in the rest frame of the B meson. We set niu = 0, M^ = 0, Mb = 
5.3 GeV, as = and Sp = 0.006. 

lb. The b —^ c double differential decay rate dT /dEedqo vs Ee for various 
values of go in the rest frame of the B meson. For the solid lines we use 
me = 1.5 GeV, Mb = 5.3 GeV,{Mx)min = Md = 1.86 GeV, as = and 
Ep = 0.006. The dashed lines have the same parameters except for {Mx)min = 
nie = 1.5 GeV. For go = OAMb the two curves coincide. 

2a. dT/dEe for semileptonic B meson decays in the rest frame of the B 
meson. The values for the parameters are rUu = 0,mc = 1.5 GeV, Mb = 
5.3 GeV, M^ = 0, Md = 1.86 GeV, «« = and various values of Ep as shown. 
For the b ^ c decays the low and high lines correspond to {Mx)min = Md 
and {Mx)min = rUe, respectively. 

2b. The b ^ u spectrum in the ACM model and the parton model. Both 
curves include QCD radiative corrections with as = 0.24. 

3. The b ^ u double differential decay rate dT/dEgdMx vs Ee for various 
values of Mx in the rest frame of the B meson. The values for the parameters 
are the same as in Fig. la. 
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